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ABSTRACT 

Aims. We present new extraction and identification techniques for supernova (SN) spectra developed within the Supernova Legacy 
Survey (SNLS) collaboration. 

Methods. The new spectral extraction method takes full advantage of photometric information from the Canada-France-Hawal tele- 
scope (CFHT) discovery and reference images by tracing the exact position of the supernova and the host signals on the spectrogram. 
When present, the host spatial profile is measured on deep multi-band reference images and is used to model the host contribution to 
the full (supernova + host) signal. The supernova is modelled as a Gaussian function of width equal to the seeing. Ax' minimisation 
provides the flux of each component in each pixel of the 2D spectrogram. For a host-supernova separation greater than ; 1 pixel, 
the two components are recovered separately and we do not use a spectral template in contrast to more standard analyses. This new 
procedure permits a clean extraction of the supernova separately from the host in about 70% of the 3rd year ESO/VLT spectra of 
the SNLS. A new supernova identification method is also proposed. It uses the SALT2 spectrophotometric template to combine the 
photometric and spectral data. A galaxy template is allowed for spectra for which a separate extraction of the supernova and the host 
was not possible. 

Results. These new techniques have been tested against more standard extraction and identification procedures. They permit a secure 
type and redshift determination in about 80% of cases. The present paper illustrates their performances on a few sample spectra. 

Key words. cosmology:observations ~ supemovae: general - methods: data analysis - techniques: spectroscopic 



1. Introduction 

Observing Type la supernovae (SNe la) for the purpose of 
constraining cosmological parameters is now a mature activity. 
Large-scale projects of detection and spectrophotometric follow 
up of hundreds of SNe la from ground-based telescopes reach 
completion today ( Astier et al.|2006 Wood-Vasey et al.|2007 



and future SNe la surveys will bring this number up to several 
thousand. These current and future surveys yield and will yield 
new samples of SNe la observations of unprecedented number 
and quality. Extracting the most possible out of these exceptional 
sets of SNe la data is a challenge for the teams involved in these 
projects. 

At redshifts greater than z ~ 0.2, the problem of spectro- 
scopically identifying supernovae is more challenging than it is 
at low redshift. At low redshift, galaxies have a larger angular 
size, so at the position of the supernova (SN) the spectrograph 
slit usually contains little contaminating host galaxy light. As 
redshift increases, apparent galaxy sizes decrease and the galaxy 



Send offprint requests to: balland@lpnhe.in2p3.fr 

* Based on observations obtained with FORSl at the Very Large 
Telescope on the Cerro Paranal, operated by the European Southern 
Observatory, Chile (ESO Large Programmes 171.A-0486 and 176.A- 
0589) 



light can be comparable to, or dominate, the supernova light at 
a given position. For observations of high-redshift supemovae, 
new techniques are therefore required to efficiently recover the 
supernova signal (see e.g., Blondin at al. 2005). 

Currently, the SNe la spectroscopic sample of the Supernova 
Legacy Survey (SNLS) represents more than 500 spectra taken 
from 8- 10m diameter telescopes, both in the Southern and 
Northern hemispheres (VLT, Gemini N and S, Keck-1 and -II) 
during large observing programmes of hundreds of hours run 
from 2003 up to the present day. The spectra obtained with such 
programmes sample a large fraction of the wavelength space (de- 
pending on the instrument used) in the visible. They have rest- 
frame phase^ comprised between -15 and -1-30 days and suffi- 
cient signal-to-noise ratios (S/N) for an unambiguous identifica- 
tion in ~ 80% of cases. 

Besides their use for cosmological purposes - that is identi- 
fication of the supernova type and determination of their redshift 
for their inclusion in a Hubble diagram, these spectra contain 
valuable information to understand the nature and the diversity 
of SNe la and test the adequacy of using them as "standardis- 
able", if not standard, candles. As an example, |Guy et aL ( 2007| l 
have implemented a number of SNLS spectra in the training 



' Throughout this paper, the phase is the restframe age of the su- 
pernova in days with respect to the B-band maximum light. 
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set of the spectrophotometric template SALT2, one of the fitters 
used to fit the Ught curves of SNe la in the 3rd year SNLS Hubble 
diagram. Moreover, comparison of the spectral properties of 
these high-redshift (< z >~ 0.5 - 0.6) supemovae with their 
low-z counterparts (from, e.g., the spectral time series currently 



collected by the SN Factory experiment, see Aldering et al. 2002 
and'Mat heson et al.|200 8 ) should help in understanding possibli 



evolutionary effects in SNe la populations (Bronder et al.|2008 



Balland et al. 2006 2007; Blon din et al.||2006t |Garavini et al.| 
2007j|Foley et al.|2008 ). Another interesting particularity of the 



SNLS spectral sample is that, due to the high average redshift 
of the SNLS SNe la, the ultraviolet (UV) part of the supernovae 
spectra down to ~ 3200 A is accessible for a large number of su- 
pernovae at various phases (see, e.g., ElUs et al. 2008 for a study 



of the UV features of 36 SNLS-Keck spectra). This is of par- 
ticular interest as it has been suggested that possible evolution 
with redshift, due to different progenitor metallicities at low and 
high redshift, might be imprinted in the spectral features of this 
region of the spectrum (Ho eflich et al.|1998 Lentz et al.|2000[ l. 
Optimal extraction of SNLS spectra and their identification are 
thus not only a crucial step in building the Hubble diagram but 
are also crucial for their subsequent use in studies of the physics 
of SNe la. 

Standard analyses of spectra are often based on the Home's 
extraction method ( |Horne|[T986 1 that uses the spatial profile of 
the source light to perform an optimal (minimum-variance un- 
biased) spectral extraction. A drawback of this approach is that 
only mild geometrical distortions due, for example, to flexure 
of the instrument are corrected for (however, see |Marsh| ( |1989| l 
for an extension of Home's method suited for highly distorted 
spectra). Another major problem of this approach is that it does 
not permit the separate extraction of the supernova from the host 
signal, except if the two components (SN and host galaxy) are 
well separated on the spectrogram (and in this case, correlations 
between the host and the SN spectra are usually unavoidable, 
even if the extraction window is adjusted to maximise the su- 
pernova signal and minimise the host contamination). In most 
cases, the SN/host separation is performed a posteriori by fitting 
a two-component model (built from a supemova and a galaxy 
template), the contribution of each com ponent to the total model 
being evaluated by & minimisation ( Howefl & Wang||2002 



IHowell et al. | |2005i [BaUand et al.1|2006| |2007| l. As the mode 
is built from a (limited) set of templates, the SN and host spec- 
tra obtained by this procedure are strongly model-dependent and 
provide at best a hint of the host contribution to the full spectral 
model. Ellis et al. (2008) have developed an improved separation 
technique based on fitting synthetic galaxy spectra to the host 
photometry measured on reference images. By doing so, they 
still rely on a spectral model to estimate the host contribution to 
the full signal. [Blondin et ar] (2005 1 use a convolution technique 
to simultaneously recover a host-free PSF-like SN component 
and its background galaxy spectrum. This method does not use 
any spectral template for the host modeling and thus improves 
over more classical extraction techniques. However, this tech- 
nique is sensitive to the width of the spatial resolution Gaussian 
kernel used to recover the background (host) component. The 
width has to be tuned by the user and highly depends on the spa- 
tial extension of the host galaxy. Moreover, the method does not 
use photometric information that potentially goes along with the 
spectrum. 

To circumvent these drawbacks and improve over standard 
extraction methods, we have developed a dedicated reduction 
and extraction pipeline, called PHASE (PHotometry Assisted 
Spectral Extraction). This new technique uses photometric pri- 



ors from Canada-France-Haw ai' Telescope (CFHT) detection im- 
ages obtained with megacam ( Boulade et al.|200? i to derive the 
exact trace of the supernova on the spectrogram. This trace is 
then used as a guide for extraction. Here, "true" photometric in- 
formation of the host contribution in each pixel, obtained from 
CFHT Legacy Survey (CFHTLS) deep reference images in var- 
ious photometric bands, is used, when possible, at the stage of 
the extraction. In the favourable cases of large enough host/SN 
separation, we do not use a spectral template to model the host 
galaxy contribution. The spatial profile of the host, measured on 
the reference image along columns parallel to the slit, is used 
instead. 

Technical choices sometimes different from the standard 
ones have been investigated to select the most efficient and ro- 
bust ones necessary for a clean and final extraction. We have au- 
tomated the reduction pipeline as much as possible to reduce the 
operator work load as well as unavoidable human errors. Besides 
the technical aspects of this new approach, in all the reduction 
and extraction processing, one of our main concerns is to control 
and propagate the noise level estimation from raw data up to the 
final reduced spectrum in the cleanest way possible. In particu- 
lar, re-sampling of data that introduce correlations among pixels 
is kept to a minimum and is done as late as possible in order 
to preserve as long as possible the independence of pixels. This 
allows us to extract a clean, optimaj^ signal-to noise spectrum, 
and in about 70% of the currently treated spectra (3rd year VLT 
spectra of the SNLS, see Balland et al.|2008l l, this is a host-free 
supemova spectrum. We discuss in this paper the improvement 
in terms of S/N over more standard extractions. 

Identification of the supemova type is then performed by 
combining spectral and light curve informations using the spec- 
trophotometric template SALT2 developed by Guy et al. ( 2007| l. 



In this approach, all available information on a given supernova 
(both photometric and spectral) is used to ease the identification. 
This contrasts with the more straightforward method of a two 
component spectral model. As underlined in, e.g.. Hook et al. 
( |20^ ; |Lidman et al.| ( |2005] l; |Balland et aT] ( |2006||2007| l, a diffi- 
culty of this approach is the possible confusion of SNe la spectra 
past maximum with Type Ic supernovae (SNe Ic) at an earlier 
phase. Using the photometric phase as a constraint is crucial to 
alleviate this degeneracy. 

In this paper, we present the techniques developed for the 
spectral extraction and the identification of our spectra. We il- 
lustrate our results with a few sample spectra obtained at VLT, 
as part of two large spectroscopic programmes running from 
June 2003 to September 2007, that are being processed using 
this technique ( Balland et al.|20()8) . In Sections 2, 3 and 4, we 
describe the PHASE technique and test it on simulated data. 
In Section 5, we illustrate PHASE results and discuss the im- 
provements over standard extractions. In Section 6, the identifi- 
cation method of SNLS spectra using the SALT2 spectral tem- 
plate is described and results are discussed on a few examples. 
Discussion and conclusion are in Section 7 and Section 8, re- 
spectively. 



2. 2D spectrum calibration using PIHASE 

Spectral extraction usually requires human interaction to define 
the locus on the frame of the object to be extracted/measured. 



- In fact, our extraction is not strictly optimal because of the PSF 
modelling as a Gaussian (see below). However, we are close to opti- 
mality (within a few percent of minimum variance). In the following, 
the word 'optimal' stands for 'close to optimal'. 
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Spectra are often distorted along the dispersion axis, and the ex- 
traction algorithm must be sophisticated enough to find and fol- 
low the source trace along this axis. Technically, one first speci- 
fies the source location as a set of pixel coordinates that is sub- 
sequently coiTected during the extraction process. 

Such a procedure, if relevant in the case of resolved sources, 
is not wanted for overlapping sources such as a point-like super- 
nova on its extended host galaxy. Beside the fact that both the 
supernova and its underlying host light are extracted together , 
the position and slope of the supernova trace itself can not be 
easily retrieved. Indeed, usual algorithms such as |Horne ( 1986| l 
seek the flux centroid, and might end up pushing the position of 
the supernova towards the galaxy brighter core. 

In order to avoid these pitfalls, we use the fact that SNLS 
produces deep images of the regions in which supernovae ex- 
plode. The supernova location and magnitudes in four photomet- 
ric bands (gM,t'M,iM,ZM) are known. Consequently, given this 
photometric information, it is possible to locate the supernova 
on the spectrogram at sub-pixel accuracy, to trace its spectrum, 
and to perform an optimised extraction able to separately mea- 
sure the supernova and host galaxy fluxes along the CCD frame. 

2.1. FORS1 long slit spectra calibration 

During the first VLT large programme (2003-2005), SNLS spec- 
tra have been acquired in long slit mode (LSS) with the FORSl 
instrument in service mode. A preliminary assessment of the 
spectroscopic data was performed within a few days in order 
to confirm the validity of the candidate. For this "real-time" 
analysis of the spectra ( Basa et al.|2008[ l, the SNLS team used 
a combination of a dedicated reduction pipeline based on the 
ESO-MIDAS software, a proper Home ( 1986 ) extraction pack- 
age and a template supernova and galaxy fitting software (SN- 



fit, Sainton 2004 



Balland et al. 2006 1 to perform an on-line iden- 



tification and redshift measurement. 

The FORSl detector is a 2kx2k CCD. To acquire supernova 
spectra, we used it with the standard collimator and mainly with 
the 300V grisrrj^ along with the GG435 order-sorting filter. The 
pixel scale is 0.2" along the spatial axis (F), and 2.65 A along 
the dispersion axis (X). PHASE has initially been implemented 
to treat the long slit spectra taken on FORS 1 . The multi-object 
(MOS) mode is cuiTently marginally supported (the majority of 
standard stars were observed in MOS mode) but its full support 
is in progress and will be used to extract the MOS spectra, mostly 
obtained during the second SNLS-VLT large programme (2005- 
2007). 

In the following, we focus on the main features of the data 
processing. Technical details can be found in Baumont (2007). 

For calibration purposes, we built bias master frames as me- 
dian of bias frames grouped per trimester (~ 50 frames). The 
scan subtracted r.m.s. of each pixel value was clipped from 8cr 
deviant events, where cr is the read-out-noise (RON) level taken 
from the raw frame header. With this threshold, most cosmic rays 
were efficiently removed. 

Normalised flat-fields were produced by grouping frames per 
run (per lunation). After an iterative median with a 5cr clip- 
ping suitable for removing most of bad pixel^ the internal lamp 
spectrum was estimated in 16 spatial bands along the Y axis and 
used to produce normalised flat-fields that reflect sensitivity fluc- 
tuation patterns in those regions. 

' And sometimes the 3001 grism for the farthest supernovae, z > 0.8 
This time, cr also includes the photon noise 



The dispersion function for science calibration was com- 
puted using a least square technique. The pixel position Xjj of 
14 isolated emission lines of Hg, He and Ar was measured in 16 
regions and averaged along the Y axis. Those lines were selected 
as the brightest and least blended ones. A set of nine dispersion 
coefficients (4th order in X', 2nd order in Y' and in XY') were ad- 
justed to the Xij values. This set turned out to be adequate for re- 
producing the dispersion solution. The 4th order coefficient Cj4 
along X has been added (it is not used by the Real-Time ESO 
pipeline) to further reduce the residuals. The residual r.m.s. is 
typically 0.2 A and reaches 0.8 A at the borders of the frame. 

We chose to use a single response function per UT (FORS 1 
moved from UTl to UT2 on June 2005) by combining all indi- 
vidual response functions derived for every standard star obser- 
vation. With this procedure, the night to night sky transparency 
variation is not corrected for. We prefer to use a uniform, easy 
to control, average calibration for the full set, rather than cali- 
brate an observation using the response function of a different 
night (ESO does not perform standard star observations for ev- 
ery night). As SNLS spectra are primarily used for the purpose of 
identification, we do not need absolute spectrophotometric cali- 
bration. Second order residual flux was sometimes found for the 
bluest standard star spectra. In this case, we estimated it and ac- 
cordingly corrected the response function. The average H2O and 
O2 atmospheric absorption spectra were computed in appropri- 
ate bands ([5800:6000], [7100:7350], [8100:8400], [8900:9850] 
A for H2O, and [6250:6350], [6800:7000], [7550:7750] A for 
O2) and used for further correction (see Section [5T| . As for the 
instrumental response, we used an average spectrum of the ab- 
sorption, computed as the average ratio of the observed spectrum 
to the reference spectrum of the standard star, multiplied by the 
instrumental response interpolated in a given absorption band. 



2.2. Combined 2D spectrogram 



In general, we took Ni„,g-2 to 5 exposures of 900s of our targets, 
with small offsets along the slit. We combined these Ni,„g frames 
of a given target to produce a clean, cosmic free sky subtracted 
2D spectrogram to be used for extraction. 

For cosmic ray removal prior to extraction, we used a tem- 
poral median filtering. The sky spectrum was obtained in 16 spa- 
tial regions by an iterative median projection along the Y axis. 
These regions are of equal width and cover the full extension of 
the CCD along the Y axis. Bright sources were detected as rows 
containing more than 5% of 5cr deviant pixels and were avoided. 
The 16 sky spectra were then averaged in wavelength space to 
yield the sky spectrum model for sky subtraction. Note that bias 
subtraction and flat-fielding were applied simultaneously during 
the first projection along Y in the iterative process. 

For a given SN candidate, bias subtraction, flat fielding and 
sky subtraction was then performed on each set of A^,,,,^ pixel val- 
ues. Then a robust, seeing-weighted average of each calibrated 
set of pixels was applied to obtain one single output pixel value. 
The seeing was computed from the active optic wavefront sensor 
system on the guiding star, recorded in the ESO logs of obser- 
vation. In rare cases (for a few standard star observations), those 
guiding star logs are lacking and the Differential IMage Monitor 
(DIMM) value was used instead. This latter value is accurate 
within 20%. 
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Figure H shows the PHASE output 2D frames for 
SN 04DlddTtaketi as an illustration of the results produced by 
the reduction technique presented above. SN 04Dldc is a Type la 
supernova (SN la) at maximum light at z = 0.21 1 well separated 
from its host core {dsN-imsi centre - 1.18")- The top panel shows 
the well separated SN (top) and host (bottom) spectra. Sky sub- 
traction residuals appear as vertical narrow bands slightly visible 
at the locus of the most luminous sky emission. Atmospheric ab- 
sorption features affecting both spectra are clearly seen at several 
places, such as the strong O2 absorption at 7600 A. The emis- 
sion and absorption features of the SN spectrum are seen as a 
modulation in the intensity of the top spectrum. The deep Si n 
absorption at ~ 6150 A is visible to the blue part of the O2 ab- 
sorption (see labels in top panel of Fig.[T]i. In the host spectrum, 
H a emission is seen as a bright spot longward of €2- Confusion 
of an emission line with a cosmic is very unlikely (a cosmic hit 
could mimic an emission line and extend spatially but it should 
be faint enough for it not to have been rejected, i.e. < 5cr of the 
average sky background). However, at high redshift (z > 0.8), 
the host flux is weak and it is more difficult to tell from the sole 
2D spectrum whether a given spot is an emission line, a statis- 
tical fluctuation, or a cosmic residual. In the present case, we 
can check in the 2D noise map shown in the lower panel of Fig. 
[T] that the observed bright spot in the spectrum is not a cosmic 
subtraction residual. Indeed, if a cosmic hit the CCD at this lo- 
cation, it would be present in this 2D noise map, which is not the 
case. This illustrates how inspection of the 2D noise map helps 
at discriminating true host lines from cosmic hits. 

3. A model for PHASE extraction 

The basic idea for the extraction from the clean 2D spectrogram 
is to define beforehand the distribution of the sources that we 
expect to find in the frame and to recover their exact position (or 
trace) on the spectrogram. To achieve this, we rely on the com- 
parison between the source spatial profiles obtained from deep 
CFHTLS reference images in three megacam photometric bands 
igM, fM and iM with Grism 300V|^ along the slit direction (the 
photometric profiles), with the corresponding profiles obtained 
from the VLT Long Slit spectrograms (the spectroscopic pro- 
files). Then, given the expected spatial profile of the sources, a 
least square linear inversion of the pixel counts on each column, 
corresponding to a given A, gives the flux of every component 
at this wavelength. In the following, we describe how we imple- 
ment this procedure in PHASE. 

3.1. Photometric priors 

To built the photometric profiles, we project the region of the sky 
covered by the slit at its position angle from the deep reference 
images in all of the megacam photometric filters covered by the 
FORSl order-sorting filter used (i.e, gM, and \m for GG435 
used with the 300V grism). This is done using SWarp (Bertin 
|& T issier'2007). This yields three ID spatial profiles that corre- 
spond to the expected galactic profiles in each filter. 

The supernova coordinates, accurately measured when 
building its light curves, are used to place the supernova com- 
ponent on these profiles. Given that supernovae are point-like 
sources, the supernova component is modelled as a Gaussian 



^ The SNLS naming scheme is the following: each supernova name 
is composed by the year of its detection (03-08), followed by the field 
name (D1-D4) and two letters in the order of discovery 

* And tm, im and zm with Grism 3001. 



profile whose width is the average seeing of the spectroscopic 
observation. As the supernova flux is not accounted for in the 
CFHTLS deep reference images, we add it "by hand" to the pho- 
tometric profiles. To do this, the flux in each band is interpolated 
from the supernova light curves at the time of spectroscopy. 

The photometric profiles are then compared to the spatial 
profiles averaged on the 2D VLT spectrogram in wavelength 
ranges corresponding to the photometric bands. For each band, 
the spatial shift between the two derived profiles (photometric 
and spectroscopic) is then computed at the maximum of their 
coiTelation function. Figure |2] illustrates the result of matching 
the VLT (solid lines) and the CFHT (dashed lines) profiles in 
gM, rM and Im (from bottom to top) for SN 04D4it, chosen as a 
test case. Correlating the two profiles yields one shift for each 
photometric band whose mean A corresponds to a given column 
^Band of the spcctrogram. Performing a linear fit of these shifts 
as a function of X yields two parameter^ 

- 1) the spatial shift between the supernova location and the 
central row of the spectrogram, at the central column, and 

- 2) the slope of the supernova location along the dispersion 
axis X. 

The maximum tolerated spatial shift is 20 pixels. In average, 
the slope is of 10"^ pixel/pixel, corresponding to a 2 pixel drift 
from border to border of the frame. 

This procedure hence allows us to measure precisely (at sub- 
pixel accuracy) the position and slope of the sources spectral 
trace on the 2D spectrogram. We now use this information to 
perform the extraction. 

3.2. Building a multi-component model 

For each column flux Fi{Y) = F{Xi, Y) of the 2D spectrogram, 
we build a A^-component spatial model as the sum of the pro- 
files C/t(Y) of individual objects in the slit, shifted by a quantity 
<5Y, according to the previously derived trace equation. Their to- 
tal flux is normalised to unity. A flux is assigned to the k''' 
component of the model and a uniform "background" compo- 
nent is added to account for possible sky subtraction residuals. 
The model for each column / reads: 



M(Y) = ^/iC,(Y-5Y,). 



(1) 



The flux fit of each component is estimated by a minimi- 
sation where: 



X 



-n 



F,(Yj)-Mi(Yj)\2 



)■ 



(2) 



Here, crj is the error associated with the column pixel data. The 
minimisation condition is equivalent to solving the following 
equation for each column: 



2; 



Co(Y')- 



2. 



Co(Y'.)C«(Y'.) 



2./ 



C«(Y')Co(Y') 



2.- 



Ca,(Y')= 



' Note it would be possible to compute the trace curvature as a second 
order term in the fitting function, but the deviation from linearity is only 
a fraction (~ 0.1 pix) of a pixel. We prefer to over-constrain the fit since 
a given band (e.g., g^ at high redshift) might contain too little flux to 
get an accurate shift value. 
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2. 



Co(y;)f,(Yj) 



CN(r.)Fi(Yj) 



where 



Yj - 6Yi corresponds to the profiles shifted ac- 
cording to the trace inchnation. 

Much care must be taken to model the source profiles, since 
their accuracy will determine the quality of the extracted spectra. 
In particular, overlapping sources is a source of anti-correlated 
noise among components. 



3.2.1. Setting the galaxy profiles 

If the extraction window is narrow enough to ensure that no field 
galaxy other than the host is present, the galaxy profile and the 
supernova Gaussian can be used to perform an extraction with 
only two components (plus a uniform background component). 

However, using a wide extraction window is preferable for 
a better sky residual estimation. This often implies that other 
object spectra are present in the window. Even in the case of 
a narrow extraction window, close galaxy pairs are sometimes 
present. In all these cases, we need to define a multi (greater than 
2) component model. Moreover, for resolved spiral galaxies, the 
spiral arm spectrum is not the same as the core spectrum. In par- 
ticular, nebular emission lines are much stronger in the arms due 
to the enhanced star burst activity. The H ii regions, ionised by 
the strong UV radiation from young and massive stars, produce 
these nebular emission lines ( |Osterbrock| 1 989| l . Improper mod- 
eling of the arm and core structure as a single component leads 
PHASE to fail to assign the correct flux to the SN and to the 
galaxy. At the wavelengths of these nebular emissions, the host 
spatial profile does not correspond to the host profile averaged 
over the whole spectral range. In the case of an SN located in 
one arm, using the average profile to model the host would lead 
to underestimate the host flux at these wavelengths, and conse- 
quently, to overestimate the supernova flux. In such cases, it is 
necessary to model the arms and the core as two different com- 
ponents. Considering that high-redshift host galaxies are usu- 
ally unresolved above a redshift of 0.7, we distinguish 3 types of 
galaxy profiles: 

- PSF : unresolved, point-like galaxies. 

- EXT : extended, regularly shaped profiles, e.g. ellipticals. 

- Mix : irregularly shaped profiles, e.g. galaxies with spiral 
arms. In practise, this profile is modelled as a PSF plus an 
EXT profiles. 

In PHASE, cuts on flux, galactic compactness, exten- 
sion minimum level and colour variation between the central 
Gaussian core (at seeing width) and the possible extension (e.g., 
arms) can be adjusted to favour one type over the others and thus 
adapt to the specific extraction case encountered. Nevertheless, 
default cuts are set to permit an automatic treatment of most 
spectra. Their values are chosen to yield the host type model 
(Mix, EXT or PSF) giving the best extraction (avoiding host 
contamination and degeneracy between the SN and host compo- 
nents) in most cases. For extended sources, we use the "bolomet- 
ric" profile (that is the sum of the galactic profiles in all covered 
filters) as a model. In order to use a noiseless extraction profile, 
it is filtered with a Gaussian of width half the seeing, oversam- 
pled at a tenth of a pixel to ease the subsequent shifting of the 
profile to follow the trace slope. Note that using the projected 
ID profile is a source of confusion when close pairs of galaxies 



happen to connect once projected along the given polar angle. 
However, sources are well defined for the subsequent extraction 
in ~95% of the case^ and this is not a major concern. 

3.2.2. Modelling the SN component 



As mentioned in Section 3.1 



the candidate supernova is an ad- 
ditional point-like source at the estimated position and of width 
equal to the seeing. The seeing variation with wavelength was 
estimated from standard star observations, as a power law of in- 
dex -0.3. This corresponds to a ~ 20 % variation of the Full 
Width at Half Maximum (FWHM) from 4000 A to 9200 A. Note 
that our power law value differs from -0.2, the value expected 
when assuming Kolmogorov turbulence (Schroeder 19871 but 
is in agreement with the value found by 'Blondin et al. ( 2005| l 
on FORS 1 spectra. Note that we use this law for a distant PSF 
host as well. However, applying the same correction to extended 
sources would require a reference image of better seeing, and 
would only have a minor effect on the profile, so we do not per- 
form it. 

To avoid degeneracy between two overlapping point-like 
sources (e.g., the supernova and an unresolved host galaxy), a 
cut of a fifth of the seeing is imposed on the separation of the 
two point-like sources. For separation values lower than this cut, 
the SN component is merged to the unresolved galaxy and no 
separate extraction is possible. As expected, this situation occurs 
with increasing frequency at higher redshifts and only the sum 
of the supernova and of the host spectra can then be obtained. 
This situation occurs in about 30% of cases (see [Balland et al.| 
2008 ) and is further referred to as the "SNGAL" case. However, 
the faintest hosts are not identified as sources since their level 
does not reach the detection threshold fixed in PHASE, and their 
spectra will only slightly contaminate the SN signal. When such 
a case occurs (or when a clear separation between the SN and 
the galaxy exists), we refer to it as a "SN"' case. 



4. Testing the validity of PIHASE extractions 

In order to assess the validity of PHASE extractions, we have 
performed a series of simulations to test the impact of vari- 
ous parameters (supernova/host separation, phase and redshift 
of the supernova, seeing) on the recovered flux for each compo- 
nent. We start from |Nobili et al. (2003 1 SNe la template spectra 
with phases ranging from -10 to H-10 days relative to maximum 
light. We assume that the supernova host is resolved and use 
the Kinney et al.| (!T996) Sb template to model its spectral en- 
ergy distribution. We follow [Blondin et al. ( 2005 i and generate 
the input synthetic 2D spectra (combined and sky background 
subtracted, for 3 exposures of 750 seconds) within the FORS 1 
instrumental environment. The profile of the host galaxy on the 
2D synthetic spectrogram is fixed from the angular size of NGC 



6181 (Kennicutt 1992) and using, as in Blondin et al. 



the [Ratnatunga et aL \ 1999 1 surface brightness profile. The flux 



2005 1, 



of the host template is scaled to the one of NGC 6181 at a dis- 
tance of 34 Mpcs. A model for the sky background (taken from 
a real FORSl exposure) is used to model the noise of the in- 
put 2D spectra (in addition to the RON and to the shot noise of 
the sources). A curvature of the spectral trace (not accounted for 
during PHASE extraction) is modelled, introducing a 0. 1 pixel 
offset from centre to border of the dispersion axis, consistent 
with the FORSl curvature measured on standard star exposures. 



** We refer here to the whole set of VLT SNe la spectra extracted 
using PHASE and presented in Balland et al. ( 2008 1 
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Figure [3] shows the input host and sky spectra, as well as the 
supernova spectrum at maximum light, for z - 0.5. 

We extract the signal with PHASE, imposing the photo- 
metric model. We perform 3 sets of simulations varying either 
the redshift (0.2 < z < 1), the SN phase (-10 < < +10 
days) or the seeing (0.4" < FWHM < 1.4"). In each case, the 
host/SN separation dsN/Host centre varies between and 2". The 
non-varying parameters are set to the default values z° - 0.5, 
d^Lim , , = 1" and FWHM'^ = 0.75". 

S N/Host centre 

To evaluate the quality of the flux recovery, we adopt the 
criterion of Blondin et al.| ( |2005| l based on the residual 6F - 
Wouf-FinM ^iF,o, i„), where Fo„, is the recovered supernova flux. 
Fin the input supernova flux, and F,„, ,„ takes into account the su- 
pernova+underlying host+sky background. If SF < 1, the input 
signal has been restored to the statistical noise limit. Figures [4]to 
|6]present 5F as a function of redshift, phase and seeing. We see 
in Fig. |4] and |5] that PHASE recovers the correct flux {5F < 1) 
for supernova-host separation dsN/Hosr centre greater than 0.5", 
even for high redshift (z > 0.9) and for early or late phases. 
For z > 0.5 and dsN/Host centre < 0.4", the galaxy becomes un- 
resolved, so that for low separation the degeneracy between the 
two component profiles is high. In fact, in such cases, PHASE 
would extract the two components together ("SNGAL" case, see 
Section [3.2.2| i. 

Note that the effect of the SN phase is very marginal (Fig. 
|5]l , indicating that low contrast of the SN relative to the host 
is not critical to recover the SN spectrum. From inspection of 
the residuals, we find that the unaccounted curvature is the main 
source of contamination of the supernova spectrum by the host 
spectrum, which will scale with this contrast. 

From Fig. l6] we see that for dsN/Hosi centre under 0.5" and a 
seeing larger than 0.9", 6F is above 1, but remains under 1.15 . 

We conclude that in a wide range of redshifts, seeing and 
phases, PHASE efficiently recovers the fluxes of the SNLS VLT 
supernova spectra. At z = 0.5 and for our median seeing of 0.7", 
PHASE is able to extract separately the supernova from its host 
provided their separation is greater than ~ 0.2" (1 pixel). 

5. Results 

5. 1 . Examples of PHASE extractions 

In Fig.|7] we present two examples of PHASE models. 

The left panel shows the "bolometric" profiles from the 
CFHT deep reference images (black dashed curve) and from the 
VLT LSS spectrum (red soUd curve) of SN 03D4ag, an SN la at 
z - 0.285 spectroscopically observed ^ 9 days before maximum 
light. The host is modelled as a Mix model (see Section [3.2. l| i 
divided into a Gaussian core (represented by a blue arrow) and 
an extended asymmetric arm profile (blue solid line). The su- 
pernova is represented as a magenta Gaussian. We thus use a 
3-component model (plus a uniform background component) to 
extract the spectrum of SN 03D4ag. 

The right panel of Fig.|7 shows the same "bolometric" model 
in the simpler case of SN 04Dldc. Here, the host is an extended 
early spiral galaxy with no fine structure, in contrast to the host 
of SN 03D4ag. It is simply modelled by the profile measured on 
the CFHT deep reference image. 

In both cases, a separate extraction of the supernova and 
the host is possible as their separation is above the default cut 
of 1 pixel ( ds N/ Host centre=^ -09" for SN 03D4ag and 1.18" for 
SN04Dldc). 

After the spectra are extracted, the flux is corrected for ex- 
posure time, instrument response, atmospheric extinction - in- 



cluding absorption from molecular oxygen and water vapour We 
also correct for slit losses under the assumption that the super- 
nova is centred in the slit and that the PSF along the dispersion 
axis has the same properties as along the spatial axis (a Gaussian 
of FWHM equal to the eff'ective seeing, with a wavelength vari- 
ation following a power law of index -0.3). Note that no interpo- 
lation of contiguous pixel values is made, neither to filter cosmic 
hits, nor to re-sample the frames at uniform wavelength bins. 
This allows an accurate estimation of the noise level, and does 
not add any further correlation among pixels. One consequence 
of this procedure is that the extracted spectra have non-uniform 
wavelength bins. 

The top panel of Fig. [8] presents the extracted spectrum of 
SN 04Dldc (red line) along with its host galaxy (blue line). 
Spectra have been rebinned for visual convenience. The host 
spectrum exhibits the features of a typical early spiral (Sa/Sb) 
with [On], CaH&K, Ha, Nii and Sii. Atmospheric absorption 
subtraction residuals are clearly seen, both in the host and the 
supernova spectra around 7600 A. 

5.2. PHASE performances and limitations 

In order to estimate its performances, we have applied PHASE 
to the extraction of all spectra of the ESO-VLT 1st large pro- 
gramme of the SNLS. The full results for the confirmed SNe la 
will be presented in |Balland et al.| ( |2008| ). A total of 829 LSS sci- 
ence and 2169 (896 MOS and 1273 LSS) calibration exposures 
have been processed to produce 208 calibrated spectrograms of 
SN candidates. In about 90% of cases, the default cuts yield an 
optimal extraction of the signal (including cases for which a sep- 
arate extraction of the SN and host components was not pos- 
sible), improved over standard extractions using classical tech- 
niques. The remaining 10% of cases are events for which ex- 
traction cuts are not suitable and thus the component definition 
is not adequate. These cases are tagged after inspection of the 
output spectrum and the residual spectrogram obtained in a first 
automatic extraction with default cuts. We treat these events on 
a case-by-case basis, adjusting the cuts to the specificity of the 
situation encountered. 

PHASE performs the extraction of all the objects present in 
the slit. The supernova spectrum is extracted separately from the 
host whenever possible (about 70% of cases). Comparison to 
standard extractions (see Section |5.3[ ) illustrates the efficiency 
of PHASE extraction within the limits of validity of the assump- 
tions made. These latter include: 

- that the coordinates of the supernova are accurate. 

- that no flux from the supernova is present in the reference 
images. 

- that the PSF is a Gaussian of width equal to the seeing. 

- that the reference image seeing is equal or close enough to 
the combined spectrogram seeing. 

Any deviation from these assumptions results in weak flux 
losses, higher noise and supernova spectrum contamination by 
its host. In particular, variation of the PSF with wavelength 
(other than the seeing dependence that we take into account) 
might be a concern. We have checked that our simple modelling 
of the PSF as a Gaussian profile yields an extraction within only 
a few percent of minimum variance. Thus, PHASE extractions 
are close to optimal. 

In PHASE, the profile shape and centroid must match as 
closely as possible those of the spectrogram. This means that one 
needs a deep reference image of seeing similar to the spectro- 
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scopic observation. However, the average seeing of the deep ref- 
erence images is 1", while spectroscopic observations at Paranal 
have an average seeing of 0.8". In this case, the exact shape of 
the source is not well reproduced by the Gaussian model (the 
galaxy is more peaked in the spectrogram than on the refer- 
ence image). Position and profile inaccuracies are reflected in 
the residual spectrogram obtained when the extracted fluxes of 
all components are subtracted from the combined spectrogram. 
If the seeing of the reference image is higher than the one of the 
VLT spectrogram, the residuals are positive at the centre of the 
spectrogram and negative on the edges. Nevertheless, for a usual 
0.1 pixel position accuracy and a 0.1" seeing accuracy (of aver- 
age 0.8"), the flux loss amounts to 5%, with a noise increase of 
about 25%. 

The bottom panel of Fig. |8] illustrates a situation for which 
the supernova coordinates are not accurate. Here, the Gaussian 
model for the supernova is located ~ 1 pixel too close to its host 
centre. The matching of the photometric profile to the spectro- 
scopic profile yields too much light in the host model and not 
enough in the supernova, as can be seen in the right panel of 
Fig. |7] The extraction yields negative residuals at the centre and 
positive residuals on the edges. 

In order to circumvent these limitations, the residuals can 
be used to iteratively adjust the position, slope and width of the 
SN component. The other (galactic) components can either stay 
fixed, or move together with the SN. Slope refinement applies 
to all sources and width refinement to point-like sources only. 
Iterating in this way yields an optimal extraction in a few passes. 
Figure [9] shows the improvement obtained after iterating in this 
way, in the case of SN 04Dldc. The black thin line is the aver- 
age profile of the residual shown in the bottom panel of Fig. [8] 
The thick blue line is the average profile of the residuals after 
five iterations. The host position is kept fixed. At the end of the 
iterative process, the supernova position has moved by about one 
pixel and its width has decreased by ~ 20 %. Inspection of Fig. 
[9] shows that the residual flux of the supernova has disappeared 
after iteration. Although useful, we do not systematically iterate, 
for the sake of preserving an automatic extraction. 

To improve further the residuals, one would need to refine 
the trace (i.e., use higher order polynomial coefficients) and PSF 
models (e.g., use Moffat functions for a better fit of the effects of 
seeing caused by atmospheric turbulence because of their longer 
tails). 

Finally, regarding computer time costs, processing spectra 
with PHASE is rather efficient. Creating the master flat-field 
frame is the most time expensive task of the calibration stage. 
It takes about 5 minutes of CPU time on an Intel Celeron pro- 
cessor running CERN-Linux to perform an iterative median on 
30 frames. The definition of photometric priors and the extrac- 
tion stage in itself takes only a few seconds. 

5.3. Comparison to standard extractions 

In this Section, we compare PHASE and standard extractions 
on a few examples. Standard extractions have been performed 
using the Horne ( 1986!) algorithm, in a way similar to what is 
done during the real-time processing of the SNLS-VLT spectra 
( |Basa et a i. 2008). With this method, the host and the supernova 
signals are often extracted together and the SN/host separation is 
done later when identifying the spectrum by adjusting a model 
built from template spectra. This model combines a fraction of 
supernova and a fraction of galaxy. From a purely practical point 
of view, using PHASE alleviates the user of a number of tasks 
such as the visual measure of the spatial shifts of the frames, the 



definition of the sky estimation regions, and the supernova and 
potential host extraction regions. 



5.3.1. SN 03D4dy: a supernova resolved from its host 

We compare here the spectra obtained with the standard and 
PHASE techniques, in the nearby subset of our sample, where 
the supernova is resolved from its host. This is the simplest pos- 
sible situation and we expect that the two extraction methods 
yield comparable results. Figure [TO] shows the two extractions 
for SN 03D4dy, an SN la at z = 0.6. The standard extraction is 
shown in the top panel while the PHASE extraction is presented 
in the bottom panel. The green spectra are the noise model for 
each extraction. The host is not visible in the slit and the model is 
a simple PSF Gaussian of width equal to the seeing. As expected 
in this simple case, the two extractions yield comparable results 
in the whole spectral range. To further compare the two extracted 
spectra, we define the quantity < S /N > as the signal-to-noise 
ratio per pixel, involving the statistical noise propagated dur- 
ing the processing (PHASE or standard) and averaged over all 
pixels. In the present example, we find that < S /N > is slightly 
higher for the standard extraction (< S /N >- 2.42) than for the 
PHASE extraction (< S /N >=2.26). 



5.3.2. SN 03D4ag: a supernova in tine arm of its spiral host 

We consider here the case of SN 03D4ag, an SN la that exploded 
in the arm of a spiral galaxy at z = 0.285. The definition of 
the sources is of critical importance in this case. In the standard 
procedure, two extractions are performed, one at the location of 
the supernova, the other in a region where the host signal is not 
contaminated by the SN. Then, this "pure" host spectrum is sub- 
tracted from the other to get the supernova spectrum. This sub- 
traction affects the colour of the final supernova spectrum. 

With PHASE, all sources are extracted at the same time pro- 
vided that the sources are properly modelled. Figure 1 1 shows 



the comparison between the two extractions. The standard ex- 
traction limits the size of the extraction window to prevent from 
too high host contamination and only part of the SN flux is 
extracted. On the contrary, if the host is correctly modelled, 
PHASE assigns the correct flux to the supernova (see Section 
|4] above). 

We find comparable < S /N > for the two extractions (13.1 
for PHASE and 15.5 for the standard extraction). The supernova 
spectrum extracted with PHASE is blue, consistent with its early 
phase {(p X -8 days). This is not the case for the standard spec- 
trum that looks like an SN la spectrum around maximum light. 



5.3.3. SN 04D2bt: a supernova close to its host centre 

SN 04D2bt is an SN la at z = 0.22 that exploded very close to 
its host centre (0.35" for a seeing of 0.62"). This renders the 
component separation difficult. Using a PSF model for the un- 
derlying galaxy, PHASE algorithm nevertheless succeeds in ex- 
tracting separately the supernova and the host, as appears in Fig. 
[T2] The host and supernova spectra obtained from the standard 
procedure are strongly correlated and the supernova flux is thus 
not accurate. With PHASE, both components are recovered cor- 
rectly at the expense of a higher noise level. 
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6. Using SALT2 for the purpose of identifying SNLS 
spectra 

In this Section, we turn to the identification of SNLS supernovae. 
Once an optimal spectrum of a potential supernova (sometimes 
combined with its host signal, see above) has been extracted us- 
ing the PHASE extraction method, it is essential to assess the 
type of the candidate and to determine its redshift for the Hubble 
diagram. 

6.1. Redshift determination 

The redshift determination is based on the identification of a 
few galactic absorption or emission lines (e.g., [On], [Om], 
Hydrogen Balmer series, CaH&K, 4000 A break). A Gaussian 
fit of one or more of these lines yields a typical uncertainty on 
the redshift of 6z < 0.001 over the available spectral range. 
This is similar to the redshift determination accuracy obtained 
for similar spectral resolutions in [Lidman et a l. (2005); Hook 
|etaL] ( |2005l l; |Howell et al.| ( |2005) ; iBailanJetal. | ( 12006) |2007, ). 

The separate extraction of the host signal from the SN in a 
vast majority of cases simplifies galaxy line identification and 
redshift determination. However, in about 20% of cases, no host 
signal is visible at all and redshift determination has to be per- 
formed on SN features, from a series of fits of the supernova 
spectrum using the SALT2 model ( Guy et al.|2007 see below). 
We allow the redshift to vary in a plausible range, given the su- 
pernova features, with a step of 6z = 0.005 and we determine 
the redshift from the best-fit solution. This latter determination 
yields a typical uncertainty of 6z ~ 0.01, ten times larger than for 
an identification from host lines. This is in agreement with red- 
shift determinations using correlation with a template, e.g. the 
SNID code developed by |Blondin & Tonry| ( [2007) 1. 



6.2. Identification 

The identification proceeds in two steps: 

1) We systematically perform a simultaneous fit of pho- 
tometric (light curves) and spectroscopic data, using SALT2 
with and without a galaxy template (see below). This yields a 
set of best-fit parameters whose values help at discriminating 
spectra with properties different from the average properties of 
the SALT2 SNe la training set ( |Guy et al.|2007| l. 



2) We then decide the type of the spectrum. This step is not 
automatic. We use both visual inspection of the fitted spectrum 
and the parameter values of the best-fit SALT2 model. The final 
decision is based on qualitative grounds. Nevertheless, using 
the SALT2 fit greatly helps in discriminating SNe la from other 
types, as both spectroscopic and photometric fit parameters 
are different for the two. In particular, we know the date of 
maximum, and thus the phase. Using the SALT2 fit allows us to 
alleviate the phase degeneracy between SNe la and SNe Ic (see 
below). 

In the following, we detail these two steps. 

The spectrophotometric template SALT2 ( Guy et al.||2007| l 
has been designed to fit the light curves of observed supernovae. 
As it is trained against a whole set of SNe la spectra templates, 
including both local, published SNe la templates and more dis- 
tant SNe la spectra, including 39 SNLS first year spectra both 
from the VLT ( Balland et al. 2008) and Gemini (Howell et al.| 
2005 1 telescopes, it is possible to use it for a combined fit of the 
photometric data and the spectrum of a given supernova. This 



procedure has been briefly discussed in |Guy et aL] ( |2007| l and 
is extensively used in the present work. It offers the great ad- 
vantage over a simple spectral fitting as the ones used in a lot 
of previous works concerned with SNe la spectral identification 
(Lidman et al."2005', 'Matheson et al.||2005t [Foley et al.|[2008 



Bafland et al. 2006, 20071 |Hook et al.|20051 



Howell et al.|2005 



Wo od-Vasey et al.)2 007) to explicitly incorporate, in the spectral 
identification process, information on the date of B-band maxi- 
mum light and colour This, in particular, alleviates the possibil- 
ity of confusing a Type Ib/c supernova (SN Ib/c) with an SN la, 
at least at early phases, and has yielded, in several cases, a firm 
type determination of dubious spectra. It is also a good way to 
discriminate Branch-normal SNe la (Branch et al.||1993j |2006| l 
from more peculiar objects for which SALT2 fails at providing 
a satisfying spectrophotometric model. 

The identification using SALT2 is based on a minimi- 
sation procedure from the combined fitting of photometric and 
spectral data for a given supernova. The main photometric pa- 
rameters entering the fit are the principal components xo, xi of 



the S ALT2 model (see Guy et al. 2007 [ for details) and the colour 
index c, defined as the difference between (B - V)sn and the av- 
erage < B-V > value for the whole training sample. The colour 
law of the model is obtained during the training process on the 
SALT2 training sample, along with the model components (see 
Fig. 3 of Guy et aL[[2007 for a comparison to a Cardelli et al. 
1989' law). The xi parameter is linked to the stretch of the su- 
pernova and can be interpreted as the number of standard devi- 
ations the stretch of the supernova under study is with respect 
to the whole training set. A set of spectroscopic parameters are 
fitted along with the photometric parameters and include a pos- 
sible host fraction, an overall normalisation parameter and a tilt 
parameter These two latter parameters allow us to take into ac- 
count possible errors in the flux calibration, due to, e.g., varia- 
tions from the time average response function used for flux cal- 
ibration, or differential refraction effects improperly corrected 
for. Higher order re-calibration parameters can be added at will, 
but we usually limit the possibility of spectral re-calibration to 
the two parameters described above, except for high S/N spec- 
tra, for which adding a third re-calibration parameter improves 
the overall fit by also correcting for "curvature" effect along the 
spectrum. 

Figure[T3]presents an example of the combined identification 
performed with SALT2. Recall that SN 04Dldc is an SN la at 
maximum light at z = 0.211, well separated from its host core 
(see top panel of Fig. 13 1. A separate PHASE extraction of the 
host and the supernova was possible, yielding an almost host- 
free supernova spectrum (Section 5.1 and Fig. [8|l. The bottom 
left panel of Fig. 13 shows the gM, fM, im, and zm photometric 



data, along with the best-fit SALT2 light curves models over- 
lapped. The fit is excellent in each band and the secondary max- 
ima, typical of SNe la, are clearly seen in the infrared light 
curves. The bottom right panel of Fig. 13 shows the correspond- 
ing spectral fit. The dashed red line is the SALT2 model ob- 
tained with no re-calibration, the solid red line being the same 
model but with re-calibration (using three re-calibration parame- 
ters in this specific case). In this example, re-calibration appears 
to have a weak effect on the final fit (significant only around 
4000-4500 A). This indicates that flux calibration is accurate. 
Strong discrepancies between the dashed and solid red lines ap- 
pear in some cases, hinting either towards a problem in flux cal- 
ibration, or, more probably, an inadequacy of the model to re- 
produce the spectral data. This happens in particular for spectra 
of supernovae that cannot be adequately reproduced by SALT2 
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(such as SNe Ic) as only SNe la are considered in the SALT2 
training sample. 

When it was not possible to extract separately the super- 
nova from the host signal with PHASE, the SALT2 fit incor- 
porates a galaxy template in the model ( |Guy et al.|2007[l. We use 



galaxy spectral series synthesised by PEGASE2 ( Fioc & Rocca- 
Volmerange 1997 1999| l for various Hubble types (Elliptical, 
SO, Sa, Sb, Sbc, Sc and Sd). Each series ranges from 1 Gyr 
to 13 Gyrs and offers a continuous sequence from blue to red 
within each Hubble type. We also use the spectral templates for 
Elliptical, SO, Sa, Sb and Sc types from |Kinney et al.| ( |TW6j l. The 
best-fit interpolated galaxy template and age are obtained as part 
of the minimising procedure. 

We have checked the quality of PHASE SN/Host separation. 
To do this, we have selected all the SNe la of the 1st VLT large 
programme for which a separate extraction of the SN component 
from its host was possible. We then did a SALT2 fit of these 
supernovae with a galactic component in the fitted model. For 
this specific test, we only used the Kinney et al. ( |1996 1 Hubble 
sequence (from Elliptical to Sc types). We discarded four extin- 
guished supernovae with very red colours and two supernovae 
for which no acceptable fit was obtained. We ended up with a 
sample of 68 SNe la (among which 47 were identified as certain 
SNe la, see our classification definition below). With the whole 
sample, we find that the average galaxy fraction in the best-fit 
model is 4.5 + 3.0 % with a dispersion of 25%. If we restrict 
ourselves to the sample of 47 certain SNe la, we find an aver- 
age galaxy fraction of 0.7 + 3.2 % and the dispersion slightly 
decreases to 22 %. These results confirm that residual host light 
contamination of the PHASE extracted supernovae is moderate. 

For each SNLS supernova for which both photometric data 
(in at least two bands) and a spectrum are available, we perform 
SALT2 fits using all galaxy types with two re-calibration pa- 
rameters. Light curves are fitted in the range -15 and -1-40 days, 
which turns out to be sufficient for most supernovae under inves- 
tigation. Spectra are fitted over their whole spectral range. The 
best-fitting model is for the minimum combined (photometric 
and spectral) value. Result parameters and fits are inspected 
to assess the identification of the candidate as a SN la (step 2 
of our identification process). Large xi absolute values hint pos- 
sible peculiarities (under or over-luminous SN). Large positive 
colour indexes are either due to large intrinsic or galactic red- 
denings, or possible SN Ib/c. Large negative values are found for 
Type II supernovae (SN II), for which no good spectral fits are 
obtained. Light curves and spectral fits are visually inspected. 
When strong re-calibration is necessary to obtain a good spec- 
tral fit, we track down the reason for this. Often, this is the sign 
of a non SN la spectrum. The combination of all these indica- 
tors usually leads to an optimised identification, even though the 
final decision relies on human judgement. 

Spectra of SNe la are classified according to the following 
scheme: SN la (certain SN la), SN la* (possible SN la but other 
types, in particular SN Ic, cannot be excluded given the S/N 
and phase), SN Ia_pec (peculiar SN la), SN (certain SN as the 
spectrophotometric properties match the known properties of an 
SN although the type is not clear), SN? (possible SN of unclear 
type, other objects such as variable stars or AGNs are possible), 
SN Ib/c, SN II. Note that this classification scheme slightly dif- 
fers from the classification adopted in |Astier et al. 1 (2006). The 
Hubble diagram is built only from SN la and SN la* supernovae 
or from SN la supernovae alone, underlying the necessity for a 
clean supernova identification. 



6.3. Results of the identification 

In this Section, we illustrate the performances of our SALT2- 
based identification through a series of representative exam- 
ples. The bulk of our SNe la identifications will be published in 
Balland et al. (2008). In the examples given below, we show the 
results of the spectral and light curve fits for the sake of com- 
pleteness. Spectra are presented in the observer frame and the 
restframe wavelength is reported on top of the corresponding 
Figure. The SALT2 model output parameters are summarised 
for each example in Table [T] 

6.3.1 . SN 05D2dt: a high S/N SN la at z ~ 0.6 



In Fig. 13 we have shown the case of a host free, low-redshift 
SN la at maximum light. Figure 14 shows the spectrum of 
SN 05D2dt, a a; -2 days SN la at z = 0.574 (about the av- 
erage redshift of the SNLS sample), observed at the VLT. This 
supernova lies at the centre of its host (see top left panel of 
Fig. 14 1 and no separate extraction was possible. The bottom 
left panel of Fig. 14 shows the full, host-contaminated spectrum 
and the SALT2 best-fit model (before and after re-calibration). 
The bottom right panel shows the host-subtracted spectrum, with 
the calibrated model overlapped. In both cases, the spectrum has 
been rebinned with 15 A bins for the sake of visual convenience. 
In the bottom left panel, the model galaxy is plotted as a blue 
solid line. Contribution of the host model to the full model is 
computed in each photometric band. A "bolometric" fraction of 
host is also evaluated over the full available spectral range and 
is given in Tab. [T] In the case of SN 05D2dt, a fair host con- 
tamination (43% of host template in the fitting model) is present 
and is well modelled by an early-type template. The model is 
tilted in order to adequately reproduce the spectrum, but this re- 
calibration is only moderate (^ 10 % flux re-calibration from 
A - 4500 A to A - 8500 A). Once the host contribution has 
been subtracted, the SNe la features are very well reproduced by 
the SALT2 model (bottom right panel of Fig. 14 1 and the iden- 
tification as an SN la is unambiguous at this phase, even if no 
clear Si n can be seen around restframe 4000 A due to the S/N 
level. 



6.3.2. SN 04D4ib: an SN la buried into its host signal 

We show here SN 04D4ib, an SN la at z = 0.704 deeply con- 
taminated by its host core signal. This is a ''SNGAL" case, for 
which PHASE cannot extract separately the supernova and the 
host. Figure 15 shows that the SALT2 fit is efficient at recovering 
the supernova spectrum (bottom right panel) from the full, host 
contaminated spectrum (bottom left panel). Here, a fraction of 
71% of a Kinney et al.| ( ,199 6) SO template has been subtracted. 
Host Ca II lines subtraction residuals are clearly seen in the su- 
pernova spectrum, but the broad SN features are well reproduced 
by the fit. No Si ii at 4000 A is visible, but given the phase of 
this SN {<p = +1 day), confusion with an SN Ic is unlikely (see 
below Section 6.3.4 1 and we can secure the type as an SN la. 



6.3.3. SN 04D4dw: a distant (z ~ 1) SN la* 

SN 04D4dw at z = 1.031 is the second farthest SN la of the 
3rd year SNLS sample. This is another ''SNGAL" case for which 
PHASE is not able to efficiently perform a separate extraction. 
Figure 16 shows the full spectrum (bottom left panel) and host- 



subtracted supernova spectrum (bottom right panel). Due to its 
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high redshift, the UV part of the spectrum is visible down to rest- 
frame 2100 A and is correctly reproduced by the SALT2 fit. This 
is not an obvious result, as the SALT2 training sample is rather 
poor in UV spectra ( Guy et a"Ll|2007[ ). The redshift is obtained 
from the presence of Ca n absorption lines. The blueshifted Ca n 
around restframe 3700 A and possibly some Si ii at restframe 
4000 A are visible in the spectrum of this supernova slightly 
past maximum (0 - 2.1 days). Given the rather low S/N, it is 
identified as an SN la*. 



6.3.4. SN 04D4jv: an SN Ic supernova 

The possible confusion of an SN la with an SN Ic is a con- 
cern when identifying a supernova candidate. Indeed, the spec- 
tral structure of an SN la slightly past maximum light resembles 
the one of an SN Ic a few days before maximum (compare for 

2004, 4 days 



instance the spectra of SN 2002bo, Benetti et al 



past maximum and of SN 19941, Filippenko et al. 



19951 5 days 



before maximum). This degeneracy is amplified at larger red- 
shift, as the spectral range of the observed spectrum and the S/N 
decrease with redshift and it is often difficult to confidently rule 
out the possibility of an SN Ic. Using SALT2 for the purpose of 
identification alleviates this degeneracy, as the phase is in prin- 
ciple known from the light curve fitting. The fact that SALT2 
fails at reproducing a spectrum is in itself an indication that it is 
unlikely to be an SN la, as far as the phase and spectral cover- 
age are well described by the SALT2 training set. Moreover, the 
values of the SALT2 colour value c might also help, as SNe Ic 
are redder than SNe la at the same phase. Indeed, SNe Ic tend to 
have high colour values, c > 0.3 (even if this sole result is not 
sufficient to unambiguously classify a supernova as an SNe Ic, 
as very extinguished SNe la would have comparable c values). 
Nevertheless, it has been noted that SNe Ic and SNe la spec- 
tra show the greatest similarity about one week after maximum 
light Pook et al.][2005| [Howell et al.|2005| l. In particular, iron 
absorption features are present in SNe Ic spectra in the range 
4200 - 5000 A, that resemble the spectral shape of SNe la ~ 1-2 
weeks after maximum light. Identification of later-phase spec- 
tra thus remains subject to caution, especially given the S/N of 
SNLS spectra at z ~ 0.5. A substantial fraction of supernovae 
classified as SN la* in our sample correspond to this situation 
of a ~ 1 week past maximum SN, for which a fair SN la match 
is obtained with SALT2 but an SN Ic solution cannot be confi- 
dently ruled out. The remaining SN la* are usually due to low 
S/N. 



The bottom left panel of Fig. 17 shows the best SALT2 match 
of SN 04D4jv (light curves and spectrum), an SN Ic at z = 0.228 
around maximum. If the light curve fit is visually acceptable, 
the SALT2 colour index is however very high: c=0.96 (see Tab. 
[TJ by construction the average c of the SALT2 training sample 
is 0). SAUr2 clearly fails at reproducing the SN structure seen 
in the spectrum around 4800 A and between 6200 and 7000 A, 
even with strong re-calibration allowed (compare the dashed and 
solid lines). A good match with an SN Ic template (SN 1997dq 
at maximum light, Matheson et al.|[200T l is obtained using the 
SN-fit software (see Section |2.1[ l and is shown on the bottom 
right panel of Fig. 17 This confirms that SN 04D4jv is an SN Ic. 



6.3.5. SN 05D4ar: an SN II supernova 
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In a similar way, we identify SN 05D4ar as an SN II (Fig 
Here, both the light curve and spectrum S AUr2 fits are very poor 
Both X] and c are very discrepant from the average values of the 



training sample: xi = 5 and c = 1.29, see Tab. [T] The bottom 
right panel shows the fit obtained with SN-fit using SN II tem- 
plates. The best-fit is for SN 1987A a few days after explosion 
(-12 days with respect to maximum light, \Pun et al.|1995| ). The 
Ha P-Cygni emission is well reproduced, as are most of the no- 
ticeable features of the spectrum. Combined to the obvious flat- 
ness of the light curves, this fact confirms that SN 05D4ar is an 
SN II supernova. 

7. Discussion 

We have taken advantage of SNLS having two indepen- 
dent pipelines to cross-check the identifications obtained using 
SALT2 with the ones done with more standard procedures of 
template fitting. For a given object, an independent assessment 
of the type and the redshift was performed on each side. On one 
side, the superfit code developed by Howell et al. (2005) was run 



on PHASE extracted spectra, using all possible SN type tem- 
plates. On the other side, SALT2 was systematically run on all 
PHASE spectra. Putting our efforts together, we have thus care- 
fully double-checked the identifications of a sample of about 250 
VLT, Gemini and Keck spectra. For VLT spectra, both extraction 
with PHASE and identification with SALT2 considerably im- 
prove the type determination. For Gemini spectra, not currently 
supported in PHASE, we only test the effect of using SALT2 as 
opposed to standard template fitting techniques. SALT2 identi- 
fications are straightforward and help in many cases to secure 
the type of the SN candidate. In some cases, the template fit- 
ting technique does a better job at reproducing the candidate 
spectrum. This usually happens for rather peculiar spectra not 
well represented in the SALT2 training sample. A few super- 
novae that had not been identified previously and that had been 
consequently classified as unknown objects, have proved to be 
SNe la in the light of our refined procedure. 

Identifications based on the fitting of the SN spectrum 
against a host and an SN template library ( Balland et al.||2006 



2007) [Howell et aL][2005| [Lidman et aL][2005[ ) improve over 
the traditional eye-guided identification of "typical" SNe la fea- 
tures. However, they rely on the template library quality that 
often suffers from insufficient phase and wavelength coverage 
and/or completeness, especially for non SNe la templates. Cross- 



correlation techniques used by ESSENCE (Matheson et al. 



2005) and SDSS ( Zheng et a l. 2008), though sophisticated, suf- 
fer to some extent from the same drawbacks, as they also rely on 
the comparison to a spectral external template dataset. Blondin 
|& Tonry[ ( |200 7^ show that the SNID code is able to confidently 
distinguish an SN Ic at redshift 0.5 with phase between -5 and + 
5 days from other SN types, depending on the use of priors on 
the redshift and age (see Fig. 21 of , Blondin & Tonry| ( ,2007) for 
details). However, this method has the hmitation that it can not 
subtract host galaxy. 

Using SALT2 takes advantage of the spectrophotometric 
model of SNe la built from a large collection of spectra and 
light curves of local and distant SNe la, a fraction of which be- 
ing SNLS supernovae. One of the key advantages of SALT2 
is that the model improves as more SNe la are detected and 
identified as part of SNLS and external collaboration efforts. 
Moreover, using simultaneous photometric and spectroscopic in- 
formation narrows the parameter space available to reproduce 
the SN properties. As an example, with SALT2, the phase is well 
constrained to within a fraction of a day, which allows us to al- 
leviate the possible confusion with other types, such as SNe Ic. 
Non SNe la events are detected, as they produce parameters val- 
ues or fits that deviate from the ones of the average sample. 
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As a model is compared to the light curves and spectral data, 
the limitation due to incomplete spectral and phase sampling is 
not as strong as in the standard template fitting technique. All 
supernova candidates are treated on an equal footing and us- 
ing SALT2 limits the subjective part of the identification intro- 
duced by human judgement. If SALT2 helps at discriminating 
non SNe la candidates, the current version of the code does not 
help in identifying their type. For such cases, recourse to a stan- 
dard template fitting technique is required to secure the type. 

8. Conclusion 

We have developed new techniques for both supernova spectral 
extraction and identification. These new tools have been devel- 
oped for the purpose of making use of the deep imaging obtained 
at CFHT to obtain an homogeneous set of SNLS supernova spec- 
tra. 

Considerable eff'ort has been put into extracting, in the most 
efficient way, VLT spectra in order to get a clean set of super- 
nova spectra at redshifts ranging from z - 0.1 to z ^ 1. One 
key feature of our so-called PHASE extraction is that the ex- 
traction is guided by computing the trace of the supernova on 
the spectrogram. This is done by using deep CHTLS reference 
images, in various photometric bands, to recover the profile of 
the non-transient objects present in the slit, from which a multi- 
component model is built by adding the (point-like) supernova 
flux as a Gaussian of width equal to the seeing. Fitting the model 
to the spectrogram light yields the flux of each component in 
each pixel. As we have shown, this improves the host-supernova 
separation over more standard techniques. A second key feature 
of the PHASE extraction is that it avoids re-sampling the data 
and correlating pixels along the procedure. This ensures a clean 
estimation of the error level associated with the signal used later 
for the identification. 

We use the spectrophotometric model of SALT2 ( |Guy et aL] 
20071 1 ™ help to determine the type of supernova candi- 
dates. Taking advantage of an ever increasing training set as new 
supernovae are discovered and followed up by the SNLS and 
other collaborations, the model is able to adequately reproduce 



SN la features in a wide range of cases ( Balland et al. 2008 show 
that it is true from phases as early as -10,-15 days up to a few 
weeks after maximum light and for a wide spectral range, includ- 
ing the UV region down to 2100 A for the most distant super- 
novae). The combined fit of light curves and spectrum tightens 
the constraints on the parameters that describe the supernova. 
Non SNe la candidates are identified on a case-by-case basis 
as their spectrophotometric best-fitting parameters deviate from 
the average properties of the SNe la training sample. Uncertainty 
however remains for the most host-contaminated (host fraction 
> 70-80%) and most distant {zl I) cases. 

PHASE extractions have been tested against simulations, 
and SALT2 based identifications have been cross-checked with 
template fitting identifications performed on the VLT data using 



the technique described injHowell et al. (2005 1. 



We find that using SALT2 for the purpose of identification 
improves the reliability of the type determination over standard 
techniques. For VLT spectra, PHASE extraction combined to 
SALT2 identification takes the most possible out of the data to 
secure a clean type and redshift determination. 
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Fig. 1. PHASE 2D combined spectrograms for SN 04Dldc. Top: the SN (top) and host (bottom) spectra after reduction. The main 
noticeable features in the SN and host spectra are indicated. Bottom: the corresponding 2D noise map, including the sky model. 
Spectra range from 4200 A (left) to 9000 A (right). Note the atmospheric absorption features affecting both the SN and host spectra. 
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Fig. 2. Match of the VLT (solid lines) and the CFHT (dashed lines) spatial profiles of SN 04D4it in three photometric bands (from 
bottom to top: gM, rM, im) arbitrarily shifted with respect to zero flux for the sake of clarity. The VLT profiles are obtained by 
averaging the spectrum in each photometric band. The CFHT profiles are measured on deep reference images. The correlation 
function of the two profiles is computed for each band and the spatial shifts are derived at the maximum of the correlation functions. 
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Fig. 3. The host (blue) and sky (green) spectra used in the simulations testing the vaUdity of PHASE extractions. The host is a Sb 
type from Kinney et al.| ( 19961 and the sky background is taken from a real FORSl exposure and is used to model the noise of the 
2D input spectra. The flux of the Sb template is scaled to the one of NGC 6181 (Kennicutt 1992) at a distance of 34 Mpcs and its 
spatial profile on the synthetic 2D spectrogram is fixed from the angular size of NGC 6181 using the surface brightness profile of 
Ratnatunga et al.| ( 1999 1. The supernova spectrum is modelled from the templates of Nobili et al. ( 2003| l, with phases ranging from 
10 to 10 days. Here, the supernova template at maximum light is shown (red line) as an example. 



Extraction quality, fct. of redshift and separation 
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Redshift 



Fig. 4. Residual 6F of a PHASE extraction performed on simulated data as a function of redshift. Values of SF lower than 1 indicate 
that the flux has been recovered to the statistical noise limit. Symbols show 6F values for SN/host separations ranging from 
to 2". For z > 0.5 and dsN/Host cemre < 0.4", 5F is greater than 1 as the galaxy becomes unresolved so that, for low separation, 
the degeneracy between the two component profiles is high. In practise, in such cases, PHASE would extract the two components 
together and not separately. 



14 



Baumont, Balland et al.: PHASE extraction and SALT2 identifications 



Extraction quality, fct. of pliase and separation 



+ 0" 




0.6' ' ' ' ' 

-10 -5 5 10 

Phase (days) 

Fig. 5. Same as Fig.|4]as a function of the phase of the simulated supernova 
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Fig. 6. Same as Fig. Has a function of the seeing of the simulated VLT spectrum. 
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Bolometric EXTRACTION Profiles of 04Dldc_589 
Bolometric EXTRACTION Profiles of 03D4ag_180 6000 n , , , , , , - 




Y pixel Y pixel 

Fig. 7. PHASE extraction "bolometric" profiles for SN 03D4ag (left) and SN 04Dldc (right). The case of SN 03D4ag illustrates 
a Mix host model composed of a Gaussian PSF core (represented as a blue arrow) and extended arms. The case of SN 04Dldc 
illustrates an EXT model. In both cases, the SN PSF model is represented as a magenta Gaussian. The horizontal dotted lines 
indicate the background level arbitrarily shifted for the sake of clarity. 
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Fig. 8. Top: PSF extracted spectra for SN 04Dldc. The host spectrum (blue line) has been rebinned with 10 A while the SN (red 
line) spectrum is rebinned with 25 A. Spectra have been corrected for atmospheric absorptions and are presented in the observer 
frame. Bottom: extraction residual spectrogram. Position and profile inaccuracies are reflected in the residual spectrogram. Here, 
the supernova in the model is located ~ 1 pixel too close to the host centre. As a result, the model profile is not accurate: too much 
light is extracted at the centre and not enough on the edges. This yields negative residuals at the centre and positive residuals on the 
edges. 
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Fig. 9. Average spatial profile of the residual spectrum before (black line) and after (thick blue line) PHASE iteration. Pixel in 
the abscissa corresponds to the centre of the slit. Residual scale is in ADUs/pixel. The SN and host locations along the profile are 
indicated as vertical arrows. Only the position of the supernova component is allowed to vary during the iteration process. 
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Fig. 10. Comparison of standard (top) and PHASE (bottom) extractions for the simple case of SN 03D4dy, an SN la at z - 0.6. 
Spectra have been rebinned with 10 A (black) and are presented in the observer frame. The green spectra are the noise model for 
each extraction. The statistical noise is plotted in red. The grey lines correspond to the intense atmospheric O [i] emissions and the 
grey bands to atmospheric O2 absorptions. 
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SN 03D4ag : Standard vs PHASE extractions 
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Fig. 11. Same as Fig. 10 for SN 03D4ag, an SN la that exploded in the arm of a spiral galaxy at z = 0.285. The Mix host model 
presented in the left panel of Fig.[7]is used to recover the different components. The PHASE supernova spectrum is bluer than the 
standard spectrum, as it should be for this ~ -8 days SN la . 
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SN 04D2bt : Standard vs PHASE extractions 
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Fig. 12. Same as Fig. 10 for SN 04D2bt, an SN la at z = 0.22 that exploded very close to its host centre, which renders the component 
separation difficult. Here, a PSF model is used for the host and permits a clean restoration of the two components at the expense of 
a higher noise level. 
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Fig. 13. SALT2 fit of SN 04Dldc, an SN la spectroscopically obsei-ved at maximum light. Top: negative deep reference image of 
SN 04Dldc host. The locus of the SN explosion is indicated with a "+" sign. The position of the slit is shown as green dotted lines. 
Bottom left: restframe gM,fM,iM,ZM light curves and corresponding SALT2 fits. Bottom right: PHASE extracted spectrum with the 
best-fit SALT2 model overlapped. Dashed red line is the raw SALT2 model. Solid line shows the best-fit model after re-cahbration 
using 3 parameters. The spectrum is presented in the observer frame and the restframe wavelength is reported on top. 
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Fig. 14. SALT2 fit of SN 05D2dt. Top left: negative deep reference image of SN 05D2dt. The locus of the SN is indicated as a "+" 
sign. The slit position is shown as green dotted lines. Top right: restframe gM,fMJM,ZM light curves and coiTesponding SALT2 fits. 
Bottom left: full PHASE extracted spectrum (including host signal) with SALT2 raw (dashed line) and re-calibrated (solid line) 
models. For the latter, two re-calibration parameters have been used. Bottom right: host-subtracted spectrum with SALT2 best-fit 
model overlapped. The spectra have been rebinned to 15 A. This supernova is identified as an SN la. 
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Fig. 17. SALT2 fit of SN 04D4jv. Top left: negative deep reference image of SN 04D4jv. The locus of the SN is indicated as a "+" 
sign. The slit position is shown as green dotted lines. Top right: restframe gM,rM,iM,ZM light curves and coiTesponding SALT2 
fits. Bottom left: PHASE extracted spectrum with SALT2 raw (dashed line) and re-calibrated (solid line) models. Even with strong 
re-calibration, SALT2 fails at correctly reproducing the spectrum aroun d 6500 A. Bottom right panel shows a fit of SN 04D4jv with 
an SN Ic template (SN 1997 dq at maximum light, Matheson et al.|200T i obtained using the SN-fiX software (see text). The spectrum 
in the right bottom panel is not rebinned. This supernova is identified as an SN Ic. 
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Fig. 18. SALT2 fit of SN 05D4ar. Top left: negative deep reference image of SN 05D4ar. The locus of the SN is indicated as a "-h" 
sign. The slit position is shown as green dotted lines. Top right: restframe gM, I'm, im, Zm light curves and corresponding SALT2 fits. 
SALT2 does not correctly fits these flat light curves. Bottom left: PHASE extracted spectrum with SALT2 raw (dashed line) and 
re-calibrated (solid line) models. Even with strong re-calibration, SALT2 fails at correctly reproducing the spectrum, namely the 
Ha P-Cygni feature at 7900 A (observer frame). Bottom right panel shows a fit of SN 05D4ar with an SN II template (SN 1987A, 
(p — — \2 days, jPun et al. 1995 ) obtained using the SN-fil software (see text). The spectrum in the right bottom panel is not rebinned. 
This supernova is identified as an SN II. 
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Table 1. Results of SALT2 fits used as an help for identification 







SN 








SALT2 model 






Name 


z 


Phase (days/max) 


"SN-Host centre 






Host Model % 


host^ 


;t^/d.o.f. 


04Dldc 


0.211 


-0.4 


1.18 


-1.56 


0.01 


NoGalaxy'' 





3.24 


05D2dt 


0.574 


-1.7 


0.01 


0.01 


0.08 


E 


43 


1.13 


04D4ib 


0.704 


0.7 


0.18 


0.60 


-0.10 


SO 


71 


0.61 


04D4dw 


1.031 


2.1 


0.11 


1.35 


0.02 


Sa 


68 


1.03 


04D4jv 


0.228 


-1 


1.5 


-0.30 


0.96 


NoGalaxy'' 





1.43 


05D4ar 


0.191 


1.3 


0.7 


5.00 


1.29 


NoGalaxy-^ 





1.71 



" Distance (in arcsec) of the supernova centre to its host centre in the PHASE extraction model. 
* From the SALT2 fit of the light curves. 

^ Percentage of host signal in the best-fitting model, averaged over the whole specttal range. 
'' The best-fitting SALT2 model has no host component. 



